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ABSTRACT 


An  apparatus  was  designed,  constructed,  and  tested  to 
measure  properties  of  single  bubbles  in  a  fluid  by  use  of  a 
multi-frequency  sound  field.  The  theoretical  background  on 
the  dual  frequency  method  for  obtaining  bubble  properties 
such  as  resonance  frequency,  rectified  diffusion  thresholds 
and  rates,  rise-time  sizing,  and  damping  coefficients  are 
discussed.  Specifications  and  design  of  the  device  are 
presented.  Sizes  of  single  air  bubbles  in  water  determined 
from  dual  frequency  methods  are  compared  to  rise-time 
sizing  for  radii  from  30  to  115  Jim,  these  two  methods  are 
shown  to  agree  within  1%.  Rectified  diffusion  rates 
measured  above,  at,  and  below  threshold  for  an  initial 
bubble  radius  of  50  |hm  in  air-saturated  water  over  a  period 
of  700  seconds  were  measured  and  the  results  demonstrate  the 
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I .  INTRODUCTION 


The  propagation  of  sound  in  water  is  significantly 
affected  by  the  presence  of  bubbles.  In  order  to  be  able  to 
predict  the  propagation  of  sound  in  bubbly  media,  the  dynamics 
of  bubbles  must  be  fully  understood.  The  reason  that  bubbles 
have  such  a  strong  influence  on  propagation  is  their  ability 
to  absorb  and  scatter  sound.  The  importance  of  these 
attenuation  mechanisms  is  determined  by  the  resonance 
frequency  and  damping  coefficient  of  the  bubbles,  quantities 
which  are  related  to  the  bubble's  radius  and  frequency 
response . 

A  number  of  articles  have  appeared  in  the  literature  over 
the  past  decade  or  so  pertaining  to  these  aspects  of  the 
dynamics  of  gas  bubbles  in  water.  The  majority  of  these 
studies  are  numerical  computations  of  such  things  as:  the 
frequency  response  of  bubbles  subjected  to  large  amplitude 
sound  fields,  the  variation  of  the  instantaneous  radius  of  the 
bubble  with  time,  and  the  variation  of  the  position  o  a 
bubble  in  the  levitating  sound  field  as  a  function  of  its 
radius.  Although  there  have  been  some  experimental 
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confirmations  of  these  predictions,  it  is  fair  to  say  that 
experiment  is  well  behind  theory. 

The  most  fruitful  experimental  techniques  involve 
acoustically  levitating  a  bubble,  thereby  allowing  extended 
periods  of  observation.  Two  methods  of  observation  have  been 
proven  particularly  successful.  The  first  involves  visually 
observing  the  bubble  through  a  telescope  and  measuring,  as 
functions  of  the  amplitude  and  frequency  of  the  drive  signal, 
either  the  position  of  the  bubble  in  the  sound  field  or  its 
average  radius  via  its  rise  velocity  [Ref.  1]  .  This  technique 
is  tedious  and  places  a  great  emphasis  on  the  skill  of  the 
observer.  In  addition,  the  kinds  of  information  are  limited. 

* 

The  second  technique  employs  illuminating  a  bubble  with  a 

laser  and  measuring  the  intensity  of  the  scattered  light  [Ref.  < 

2].  This  technique  takes  the  burden  off  the  observer  at  the 

expense  of  a  more  complex  experimental  setup.  In  addition, 

the  data  must  be  inverted  by  the  use  of  Mie  scattering  theory 

to  infer  the  bubble  radius.  Nonetheless,  this  technique  has 

great  potential.  In  particular,  it  can  yield  the 

in  cantaneous  radius  of  a  bubble,  a  very  desirable  quantity. 

A  different  observation  technique  is  addressed  in  this 
thesis.  It  is  based  on  a  bubble-sizing  method,  developed  by 
Newhouse  and  Shankar  [Ref.  3 ] ,  using  the  nonlinear  mixing  of 
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two  frequencies  (dual  frequency  method)  which  can  yield  the 
same  information  that  is  available  through  visual 
observations,  but  with  much  less  effort  on  the  part  of  the 
observer.  It  can  also  provide  information  about  the  frequency 
response  of  the  bubble,  which  is  not  attainable  through  visual 
observations.  In  the  present  configuration  of  this 
observation  technique,  it  is  not  possible  to  obtain  the 
instantaneous  radius  of  the  bubble. 

The  theory  necessary  to  understand  the  detection 
technique,  as  well  as  the  measurements  performed,  are 
presented  in  the  next  chapter.  Next  the  design  of  the 
measuring  device  and  the  associated  experimental  setup  are 
discussed.  The  results  of  measurements  of  bubble  radius, 
growth  rate,  and  frequency  response  are  then  presented  and 
discussed.  The  thesis  concludes  with  a  summary. 
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II.  THEORY 


The  goal  of  this  thesis  is  to  develop  the  frequency  method 
of  investigating  the  dynamics  of  bubbles.  In  particular,  it 
is  desired  to  determine  the  size  of  a  bubble,  its  growth  rate 
when  exposed  to  a  sound  field,  and  its  damping  coefficient. 
The  first  two  quantities  can  be  determined  from  measurements 
of  the  resonance  frequency  of  the  bubble.  The  last  quantity 
can  be  determined  by  measuring  the  radiated  acoustic  pressure 
amplitude  as  a  function  of  driving  frequency. 

The  principle  behind  the  dual  frequency  method  is  given 
in  the  next  section,  followed  by  a  discussion  of  the  resonance 
frequency  of  a  bubble.  Then  a  description  is  given  of  an 
alternate  method  of  determining  the  size  of  a  bubble  by 
measuring  the  time  required  for  it  to  rise  through  a  known 
distance.  Next,  the  theory  governing  the  growth  rates  of 
bubbles  in  sound  fields,  "rectified  diffusion",  is  presented. 
The  chapter  concludes  with  a  discussion  of  the  frequency 
response  and  damping  coefficient  of  bubbles. 

A.  DUAL  FREQUENCY  METHOD 

The  majority  of  research,  both  experimental  and 
theoretical,  pertaining  to  bubble  dynamics  has  considered  the 


4 


problem  of  a  bubble  excited  by  a  monofrequency  sound  field. 
The  fundamental  equation  describing  this  phenomena  is  the 
Rayleigh-Plesset  equation  [Ref.  4:eq.  4,  p.  122]  which  is  a 
nonlinear  equation  expressing  the  radius  of  the  bubble  as  a 
function  of  time. 

p/?«+-|-p/?2  =  (/’.  +  ^-)(^)3,'-(^  +  ^-)-4M-|-  +  />)COSWi,  (eqn .  1) 
where 

p  =  density  of  the  liquid 

P  =  viscosity  of  the  liquid 

a  -  surface  tension 

R0  =  equilibrium  bubble  radius 

R  =  instantaneous  bubble  radius 

Pc  =  ambient  pressure  amplitude 

Pj =  incident  sound  pressure  amplitude 

)/  =  polytropic  exponent  [l<»J<cp/cv  (ratio  of  specific  heats)] 
=  angular  frequency  of  the  sound  field 

The  dots  indicate  time  derivative.  The  derivation  of  the 
Rayleigh-Plesset  equation  assumes  that  the  bubble  can  be 
treated  as  a  highly  compressible,  spherical  body  immersed  in 
an  incompressible  fluid.  All  terms  are  expressed  as  energy 
per  unit  volume.  The  kinetic  terms  are  displayed  on  the  left 
end  side  of  the  equation.  The  first  term  of  the  right  end 
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side  represents  contributions  from  thermal  processes.  The 
following  two  terms  are  the  effects  of  surface  tension  and 
viscous  stress  of  the  bubble  surface.  The  last  term 
represents  the  energy  resulting  from  the  driving  force. 

The  method  used  in  this  thesis  involves  two  sound  fields. 
This  dual  frequency  method  is  described  by  including  a  second 
driving  term  P4  coswjt  on  the  right-hand  side  of  the  Rayleigh- 
Plesset  equation.  Newhouse  and  Shankar  [Ref.  3  gives  an 
extensive  treatment  of  bubble  response  under  the  influence  of 
two  sound  fields  of  different  frequencies,  £1  and  £2.  fl  is 
much  higher  than  the  bubble's  resonance  frequency,  but  with 
a  wavelength  in  the  fluid  still  much  longer  than  the  diameter 
of  the  bubble.  £2  is  at  or  near  the  bubble's  resonance 
frequency.  The  bubble  is  excited  into  large  amplitude 
resonant  oscillations  by  £2  while  also  responding  to  £1, 
although  with  a  much  smaller  oscillation  amplitude. 

Because  the  behavior  of  a  bubble  undergoing  large 
amplitude  oscillations  is  nonlinear,  the  bubble  radiates  sound 
not  only  at  the  incident  frequencies  £1  and  £2,  but  also  at 
the  sum  ( £l+f2 )  and  difference  (£l-£2)  frequencies.  All  other 
possible  combinations  are  also  radiated,  but  with  much  lower 
amplitudes . 
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When  the  frequency  £2  equals  the  resonance  frequency  of 
the  bubble  £0,  the  radiated  pressure  amplitude  ( £l  +  £2 )  is 
given  by  Newhouse  and  Shankar  [Ref.  3 : eqn .  21,  p.  1475] . 

P  P 

p+  — - LJ _  (eqn.  2) 

pcolR0Sr 

where 

P1  «  acoustic  pressure  amplitude  at  the  frequency  fl  at  the 
location  of  the  bubble 

P2  =  acoustic  pressure  amplitude  at  the  frequency  £2  at  the 
location  of  the  bubble 

P4  «=  acoustic  pressure  amplitude  at  the  frequency  (f l  +  f2) 
at  the  location  of  the  receiver 

a)0  -  angular  resonance  frequency  of  bubble 

<5  =  damping  coefficient 

r  «  range  from  bubble  to  receiver 

One  advantage  of  the  dual  frequency  method  over 
monofrequency  methods  is  that  the  sum  and  difference 
frequencies  are  generated  only  by  resonant  bubbles  as  shown 
in  Figure  1.  Hence,  it  is  not  sensitive  to  larger  nonresonant 
bubbles  or  solid  contaminants  in  the  water.  The  reason  is 
that  the  generation  of  the  sum  frequency  requires  that  the 
bubble  be  undergoing  large  amplitude,  nonlinear  oscillations. 
Even  though  the  backscattered  pressure  from  the  bubble  peaks, 
dips,  and  rises  linearly  as  shown  in  curve  b  of  Figure  1,  the 
sum  and  different  frequencies  (fl+fO)  or  (fl-fO) ,  are  strong 
for  only  one  bubble  size  as  shown  in  curve  a. 
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0-1  h  ft  X  -0.2 


(?03/9euXp)  ?d 


( ^mo/agu^p)  4  <3 


Bubble  in  Water  with  JfO  =  20  kHz ,  Curve 
P*  and  Curve  b:  P. . 


B .  RESONANCE  FREQUENCY 

The  concepts  of  bubble  resonance  and  resonance  frequency 
are  crucial  to  the  dual  frequency  method.  Therefore  it  is 
worthwhile  reviewing  them.  Devin  gives  the  linear  second 
order  differential  equation  of  motion  for  a  bubble  [Ref.  5:eq. 
19,  p.  1656] . 

(  )i;  +  bv  4  (y  ~~  )  =  -  /’  exiH/ni/)  (eqn.  3) 

where 

v  =  change  in  volume  fr  :>m  the  equilibrium  bubble  volume 
b  =  dissipation  coefficient 

V  =  ratio  of  specific  heats 

Vn  =  equilibrium  bubble  volume 

P  =  forcing  pressure  amplitude 

If  the  dissipation  is  negligible,  equation  3  gives  the 
resonance  frequency  of  the  bubble  [Ref.  5:  eq.  21,  p.  1656J 

r= — ! — ( )'/2  (eqn.  4) 

Jo  Anli0  1  P  ’ 

which  is  often  called  Minnaert's  expression.  Figure  2  is  a 
graphical  representation  of  equation  4  for  air  bubbles  in 
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Figure  2 .  Equilibrium  Radius  for  Air  Bubble  in 
Water 
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water.  Surface  tension  has  been  neglected  in  this  discussion 
because  its  effects  are  unimportant  for  bubble  sizes  of 
interest  in  this  study. 

C.  SIZING  07  BUBBLES  BY  RISE  TIME 

In  order  to  verify  the  dual  frequency  technique,  an 
alternate  method  of  determining  bubble  size  is  necessary.  A 
relatively  accurate  method  of  finding  the  radius  of  a  bubble 
is  to  measure  the  time  required  for  a  bubble  to  rise  through 
a  known  distance  in  still  water  under  the  influence  of 
buoyancy . 

As  documented  in  [Ref.  6],  a  bubble  rises  at  constant 
speed  due  to  the  balance  between  buoyancy  and  drag  forces. 
Since  these  forces  depend  on  radium  differently,  the  speed  of 
a  bubble  varies  with  radius . 

The  terminal  rise  speed  of  a  bubble  is  [Ref.  6] 

(eqn.  5) 

where 

U  =  rise  speed 

a  =  bubble  radius 

g  =  acceleration  due  to  gravity 

CD  =  drag  coefficient 
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The  drag  coefficient  developed  by  Schiller  and  Nauraan  was 
used  fc-  these  computations  [Ref.  6:eq.  4,  p.  6] 

CD  =  -|p-  (I  +  0.15fl°  687)  (eqn.  6) 

where 

R.  =  Reynold's  number  =  2aU/vf 

vf  =  kinematic  shear  viscosity  of  fluid 

Figure  3  displays  the  bubble  radius  relative  to  the  time 
for  a  bubble  to  rise  25  mm  in  water.  This  method  is  not  valid 
for  bubble  radii  smaller  than  15  flm  since  the  size  is  too 
small  to  be  observed. 

D.  RECTIFIED  DIFFUSION 

In  order  to  determine  the  utility  of  the  dual  frequency 
techniques,  one  measurement  we  can  make  is  of  the  growth  rates 
of  bubble  exposed  to  a  sound  field.  The  sound  field  will 
cause  the  bubble  to  grow  (or  dissolve)  by  a  process  known  as 
rectified  diffusion.  Crum  [Ref.  7]  gives  a  detailed 
discussion  of  the  subject.  Rectified  diffusion  involves  the 
slow  growth  of  a  pulsating  gas  bubble  .due  to  an  average  mass 
flow  into  the  bubble  as  a  function  of  cime.  The  growth  of  a 
bubble  can  be  detected  by  the  associated  change  in  its 
resonance  frequency.  The  threshold  acoustic  pressure  for 
growth  of  a  gas  bubble  is  [Ref.  7:eq.  19,  p.  217) 
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where 

a)  =  angular  frequency  of  excitation 

ct  =  concentration  of  dissolved  gas  in  the  fluid  far  from 
the  bubble 

ce  =  equilibrium  concentration  of  gas  in  the  fluid 

K  «  thermal  conductivity  of  the  gas  in  the  bubble 

Pth  =  threshold  acoustic  pressure  amplitude  at  excitation 
frequency 

Table  1  displays  some  results  using  equation  7. 

Crum  [Ref.  1]  also  discusses  the  equations  that  describe 
the  change  in  size  of  an  air  bubble  that  is  present  in  a 
liquid  and  experiencing  a  sinusoidal  fluctuation  in  pressure. 
The  expression  for  the  rate  of  change  of  the  bubble  radius 
with  time  is  given  by  [Ref.  l:eq.  13,  p.  205] 
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(eqn.  9) 
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TABLE  1. 


THRESHOLD  ACOUSTIC  PRESSURE  AMPLITUDES  FOR 
RECTIFIED  DIFFUSION  OF  VARIOUS  AIR  BUBBLES  IN  AIR 
SATURATED  WATER  DRIVEN  AT  RESONANCE 


BUBBLE  RADIUS 
(  m) 

RESONANCE  FREQUENCY 
(kHz) 

THRESHOLD  PRESSURE 
(Pa) 

100 

31.4 

1440 

50 

62.1 

5690 

10 

303.8 

8680 

5 

618.0 

1.515  MPa  * 

*  carrier 

frequency  =  2.65  mHz 
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3  rnR0 


where 

D  =  diffusion  constant 

P*  =  applied  acoustic  pressure  amplitude 

Pq  =  instantaneous  pressure  amplitude  of  the  gas  within 
the  ubble 

Bubbles  can  either  grow  or  dissolve  depending  upon  a 
variety  of  conditions  as  shown  in  Figure  4.  Notice  that  it 
takes  hundreds  of  seconds  to  cause  a  major  change  in  bubble 
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Rubble  Radius  (|im) 


Figure  4 .  Rectified  Diffusion  Measurements  by  Crum, 
Solid  Lines  Representing  Theory  [Ref .  7 : fig . 
4,  p.  207]. 


radius.  Therefore,  as  long  as  the  pressure  amplitude  used  in 
the  dual  frequency  method  are  at  or  below  threshold  amplitude 
and  they  are  not  applied  for  long  periods  of  time,  the  dual 
frequency  technique  should  not  significantly  affect  the  size 
of  the  bubbles  being  investigated. 

B.  FREQUENCY  RESPONSE,  DAMPING,  AND  QUALITY  FACTOR 

The  frequency  response  and  damping  coefficient  have  been 
two  of  the  most  difficult  characteristics  of  bubbles  to 
accurately  measure. 

The  total  damping  constant  originates  from  three 
processes : 

-  thermal  damping  $  th  due  mainly  to  the  thermal  conduction 
between  the  gas  in  the  bubble  and  the  surrounding  liquid 

-  Sound  radiation  damping  ^rad 

-  Viscous  damping  d  vls  due  mainly  to  viscous  forces  at  the 
gas-liquid  interface. 

The  total  damping  constant  is  the  sum  of  the  damping  due 
to  these  three  processes: 


<5  -  <\fi  +  ^rad  +  &  vis 


(eqn.  15) 
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Prosperetti  [Ref.  8]  gives  an  elaborate  discussion  of  the 


three  terms.  They  are  [Ref.  8:eq.  29+40,  p.  72] 
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with 


Kg 

Xg  =  — p—  (eqn.  19) 


sinkO  ±  sin  0 
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(eqn.  20) 


where 

c  =  speed  of  sound  of  fluid 

P,(R„)=  pressure  acting  on  the  inner  side  of  the  bubble 
surface 

0  =  ratio  of  the  bubble  radius  to  thegas  thermal  diffusion 

length 

Xq  =  thermal  diffusivity  of  the  gas 
Kq  =  thermal  conducivity  of  the  gas 
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P9  =  density  of  the  gas 

0^  =  specific  heat  at  constant  pressure 

For  the  low  frequencies,  equation  16  shows  that  viscous 
damping  dominates  at  small  radius.  In  contrast,  equation  18 
shows  that  thermal  damping  dominates  at  larger  radius . 
Equation  17  shows  that  the  acoustic  contribution  is  negligible 
at  low  frequency.  The  reader  is  referred  to  [Ref.  8]  for  a 
more  complete  discussion. 

For  oscillations  of  a  bubble,  the  quality  of  factor  Q  is 
[Ref.  5:eq.  23,  p.  1657] 

Q  =  ~~r  (eqn.  21) 

J  2  J 1 

where  fl  and  f2  are  the  two  frequencies  below  and  above 
resonance  at  which  the  average  sound  power  radiated  by  bubble 
has  dropped  to  half  its  resonance  value.  The  total  damping 
constant  is  also  defined  as  the  reciprocal  Q  of  the  bubble, 
<5=  1/Q.  Measurements  of  the  bubble  response  can  therefore 
yield 
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III.  APPARATUS  AMD  MEASUREMENTS 


The  first  step  toward  accomplishing  the  goal  of  this 
thesis  is  to  design  an  apparatus  which  utilizes  the  dual 
frequency  method  to  study  bubbles  in  a  fluid.  The  design  and 
use  of  the  device  is  described  in  this  chapter. 

A.  DESIGN  CONSIDERATIONS 

It  is  proposed  to  use  the  dual  frequency  technique  to 
determine  the  sizes  and  responses  of  bubbles  in  the  range  10 
to  125  ftm  radius.  Since  it  is  important  that  the  diameter  of 
the  bubble  be  small  with  respect  to  the  wavelength  associated 
with  the  frequency  fl,  a  maximum  carrier  frequency  of 
approximately  3  MHz  is  necessary  to  insure  that  the  acoustic 
wavelength  is  at  least  four  times  the  size  of  the  largest 
bubble  radius  of  125  ^m.  Radii  from  125  Jim  down  to  10  Jim 
correspond  to  resonance  frequencies  from  approximately  30  kHz 
to  330  kHz.  Assuming  that  harmonics  of  fl  higher  than  the 
fifth  will  be  insignificant  leads  to  a  minimum  carrier 
frequency  of  about  six  times  the  highest  bubble  resonance  or 
approximately  2  MHz.  Therefore,  a  value  of  2.65  MHz  was 
chosen  for  fl  in  the  experiment. 
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The  frequency  of  f2  of  the  other  sound  field  will  be  swept 
from  below  to  above  the  resonance  frequency  of  the  bubble. 
The  swept  rate  must  be  slow  enough  to  insure  that  the  bubble 
reaches  its  steady  state  response  at  a  given  frequency. 
However  the  sweep  rate  must  be  high  enough  and  the  amplitude 
low  enough  so  that  the  bubble  does  not  grow  due  to  rectified 
diffusion . 

To  be  able  to  make  measurements  on  a  single  bubble  over 
a  prolonged  period  of  time,  it  is  necessary  that  the  bubble 
remain  in  the  same  place  during  the  interrogation  process. 
To  accomplish  this  requirement,  the  bubble  will  be  levitated 
in  the  sample  region.  The  levitation  field  can  also  serve  as 
the  driving  field  in  order  to  study  rectified  diffusion  and 
bubble  damping.  A  discussion  of  bubble  levitation  is  given 
in  [Ref.  4]  . 

From  the  definition  of  all  the  above  mentioned  features, 
the  use  of  the  dual  frequency  method  requires  that  the  device 
has  the  following  capabilities: 

-  A  method  of  generating  bubbles 

-  A  method  of  levitating  bubbles 

-  A  method  of  producing  a  carrier  frequency  fl  and  an 
exciter  frequency  f2 

-  A  method  of  receiving  the  pressure  signals  from  the 
bubble  and  producing  a  frequency  spectrum  of  the  received 
signal 
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1 .  Generating  Bubbles 

A  method  for  generating  bubbles  of  radii  sizes  10  to 
125  |lm  is  necessary  for  this  thesis.  A  way  to  generate 
bubbles  of  such  sizes  is  to  stimulate  rectified  diffusion  of 
microbubbles  already  present  in  the  saturated  fluid.  Once  the 
bubbles  are  at  the  desired  size,  the  amplitude  of  the  driving 
field  is  reduced,  thereby  halting  the  growth. 

2 .  Bubble  Levitation 

Since  the  position  of  the  bubble  is  important  in  this 
device,  it  is  desirable  to  hold  the  bubble  at  precisely  the 
same  position.  This  can  be  accomplished  with  the  levitation 
cell  described  below. 

From  a  vertical  plexiglass  tube  as  a  container  for 
the  fluid,  transducers  can  be  employed  to  excite  a  standing 
wave  pattern  along  the  length  of  the  tube.  A  bubble  can  be 
levitated  in  the  standing  wave  pattern  by  adjusting  the 
pressure  amplitude  to  counteract  the  buoyancy  force  of  the 
bubble.  By  selecting  the  appropriate  transverse  normal  model, 
with  maximum  velocity  (minimum  pressure)  at  the  walls  and 
minimum  velocity  (maximum  pressure)  at  the  center,  the  bubble 
can  be  held  at  the  center  of  the  cylinder.  A  frequency  that 
produces  a  normal  model  (0,  1,  4)  for  a  cylinder  (0,  r,  z)  as 
presented  in  [Ref.  9]  will  produce  such  a  result. 
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3 .  Carrier  and  Exciter 

The  two  sound  fields  required  in  the  dual  frequency 
method  can  be  generated  by  commercially  available  immersion 
transducers.  The  carrier  transducer  must  be  capable  of 
generating  sufficient  amplitudes  at  the  specified  carrier 
frequency  fl.  On  the  other  hand,  the  exciter  transducer  must 
be  capable  of  generating  signals  over  a  range  of  frequencies, 
f2,  at  constant  amplitude. 

4 .  Reception  and  Display  of  Received  Signal 

Once  the  bubble  has  been  excited,  a  transducer  is 
needed  to  record  the  resulting  pressure  radiated  by  the  bubble 
as  a  function  of  the  frequency  f2.  As  described  in  the  theory 
section,  a  bubble  will  radiate  much  higher  amplitudes  at  its 
resonance  frequency  fO  than  at  any  other  frequencies.  This 
received  signal  can  be  displayed  with  the  help  of  a  frequency 
spectrum  analyzer. 

B .  EXPERIMENTAL  SET-UP 

Following  the  establishment  of  the  design  considerations, 
it  is  now  possible  to  build  the  device.  This  section  will 
describe  the  transducers  employed,  the  actual  device,  the 
electronic  equipment,  and  the  schematic  of  the  total  system. 
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1 .  Transducers 


As  described  in  paragraph  A,  four  transducers  are 
required.  A  3-inch  diameter,  1-inch  long,  piezoelectric 
ceramic  type  is  used  to  establish  standing  waves  within  the 
plexiglass  cylinder  of  same  diameter.  Located  at  the  top  of 
the  cylinder  as  shown  in  Figure  5,  it  is  used  at  a  frequency 
of  21.90  kHz  to  excite  the  desired  normal  model  in  the 
cylinder . 

Stuck  to  the  bottom  of  the  plexiglass  cylinder,  a 
Panametrics  V301  1-inch  diameter  circular  piezoelectric 
immersion  transducer  (serial  number  93596)  is  used  as  the 
exciter.  The  selection  of  this  transducer  was  based  upon  the 
fact  that  the  frequency  response  over  the  range  of  20  to  100 
kHz  is  flat,  within  5dB.  This  transducer  could  have  been 
located  almost  anywhere  as  long  as  the  bubble  can  be  excited. 
To  demonstrate  this,  a  small  piezoceramic  hydrophone  will  be 
used  as  the  exciter  later  to  compare  results  from  two 
different  exciters. 

Two  Panametrics  V323,  0.25-inch  diameter  circular 
piezoelectric  immersion  transducers  are  used  as  the  carrier 
(serial  number  923564)  and  the  receiver ( serial  number  93566) . 
These  two  transducers  are  located  6.0  cm  from  the  top  of  the 
cylinder.  There  were  placed  in  position  by  a  plexiglass  block 
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Levitation  transducer 


Figure  5.  Transducer  Locations  for  Multi -Frequency 
Method 
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formed  to  hold  them  with  the  help  of  four  screws.  Because  the 
beamwidth  of  these  two  transducers  is  only  10  degrees  in  the 
frequency  range  used,  precise  alignment  is  necessary  to  ensure 
that  both  beam  patterns  intersect  at  the  location  of  the 
bubble . 

2 .  Apparatus 

After  consideration  of  the  requirements  placed  on 
transducer  position  discussed  in  the  last  paragraph,  the 
transducers  were  attached  to  the  cylinder  at  the  positions 
shown  in  Figure  5.  This  arrangement  enabled  the  experimenter 
to  actually  see  the  bubble,  simplifying  the  task  of  centering 
the  bubble  between  the  earlier  and  receiver  transducers. 

The  carrier  transducer  transmits  a  CW  sinusoidal 
signal  at  2.65  MHz.  To  reduce  the  amount  of  reflection  of 
this  signal,  a  castor  oil  chamber  was  placed  directly  across 
from  the  carrier  transducer.  The  castor  oil  was  kept  separate 
from  the  water  in  the  test  cell  by  means  of  a  thin  mylar 
window . 

The  brass  plate  used  to  incase  the  exciter  transducer 
also  served  as  the  bottom  of  the  cylinder  insuring  a  pressure 
antinode  at  the  bottom  of  the  chamber.  The  overall  dimensions 
of  the  actual  chamber  are  17.8  cm  high  by  7.6  cm  diameter. 
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3 .  Electronic  Equipment 


This  paragraph  will  detail  the  electronic  equipment 
shown  in  Figure  6. 

To  produce  and  levitate  the  bubble,  a  HP8904A 
Multifunction  Synthesizer  is  used  to  produce  a  21.90  kHz  sine 
wave.  This  signal  is  amplified  wiLh  an  Amplifier  Research 
model  100A15  power  amplifier  which  has  a  variable  gain 
control.  The  output  was  monitored  on  a  Tektronix  2445A 
oscilloscope . 

The  exciter  signal  is  generated  by  a  HP3314A  Function 
Generator  in  the  sweeping  model  from  20  kHz  to  100  kHz.  This 
signal  is  amplified  by  an  Amplifier  Research  model  50A15 
amplifier.  A  low  pass  RC  filter  was  used  to  eliminate 
frequency  components  above  500  kHz,  therefore  reducing  the 
changes  of  interference  with  the  carrier.  The  Tektronix  2445A 
oscilloscope  was  also  used  to  monitor  the  signal. 

The  carrier  signal  is  a  2.65  MHz  sinusoidal  waveform. 
The  signal  is  generated  by  a  HP8640R  signal  Rpnerator  which 
has  a  very  good  spectral  purity.  This  purity  is  required 
because  the  sideband  signals  product-  i  by  the  dual  frequency 
method  flank  the  carrier  on  a  spectral  plot. 

The  received  signal  is  amplified  with  a  Minicircuits 
model  ZHL-32A  preamplifier  having  a  gain  of  27dB  +/-  2dB  over 
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Figure  6.  Multi -Frequency  Method  Setup 


the  range  of  interest.  The  amplified  signal  is  then  sent  to 
a  HP3585A  Spectrum  Analyzer.  To  reduce  the  range  of  the 
required  frequency  sweep,  thereby  increasing  the  resolution, 
only  the  sum  frequency  would  be  swept  and  displayed.  For 
example,  with  a  carrier  of  2.65  MHz  and  an  exciter  sweeping 
from  20  to  100  kHz,  the  spectrum  analyzer  will  sweep  from  2.67 
to  2.75  MHz.  The  resulting  spectrum  would  then  be  plotted  on 
the  HP7090A  plotter.  The  sweep  time  of  the  analyzer  was  set 
to  match  that  of  the  exciter. 

4 .  Sample  Fluid 

The  fluid  used  in  these  measurements  is  deionized,  air 
saturated  water.  A  thermometer  is  used  to  record  the  room 
temperature  and  a  tensiometer  to  measure  the  surface  tension 
of  the  water.  Between  each  data  run,  the  water  was  vigorously 
shaken  insuring  an  air  saturated  sample. 

5 .  Calibration 

Due  to  the  complexity  of  the  acoustic  field  in  the 
cylinder  caused  by  the  three  active  transducers,  it  was  judged 
inappropriate  to  determine  the  absolute,  free  field, 
calibrations  for  each  transducer.  Introducing  a  small 
calibrated  hydrophone  and  measuring  the  actual  field  present 
at  the  location  of  the  bubble  was  a  more  feasible  and  accurate 
way  of  determining  the  acoustic  pressures  at  the  bubbles. 
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Two  hydrophones  were  used:  an  LC-10  (serial  number 
B-907)  and  a  small  transducer  built  from  a  Channel  5500  thin 
walled  cylindrical  ceramic  reversible  transducer.  A  full 
reciprocity  calibration  as  presented  in  Chapter  II  of  [Ref. 
2]  was  executed  on  both  hydrophones  producing  a  relatively 
(+/-ldB)  accurate  receiver  sensitivity  from  5  kHz  to  100  kHz. 
It  was  then  possible  to  insure  that  the  carrier  and  exciter 
would  not  produce  pressure  sufficient  to  cause  rectified 
diffusion.  When  rectified  diffusion  was  desired,  it  was 
produced  with  the  21.90  kHz  levitation  field.  A  scale  marked 
on  the  variable  gain  control  on  the  amplifiers  indicated  the 
actual  pressure  amplitudes  acting  on  the  bubble. 

6 .  Procedures 

The  bubble  radius  was  measured  using  both  the  rise 
time  method  and  the  rectified  diffusion  method.  The  former 
one  requires  that  the  carrier  and  exciter  pressure  amplitudes 
by  zero  while  allowing  the  bubble  to  rise  on  its  own  after 
turning  off  the  levitation  transducer.  The  time  for  the 
bubble  to  travel  between  two  marks  separated  vertically  by  25 
mm  placed  on  the  plexiglass  is  measured  using  a  stopwatch. 

For  the  latter  case,  rectified  diffusion  method,  the 
bubble  is  levitated  to  the  desired  position.  Next  the  desired 
pressure  amplitude  is  applied  to  the  bubble  by  the  levitation 
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sound  field.  At  approximately  every  100  seconds,  the  bubble 
radius  is  measured  using  the  dual  frequency  method  with  the 
use  of  the  carrier  and  exciter  transducers.  This  process 
takes,  on  the  average,  15  seconds. 
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IV.  RESULTS  AMD  DISCUSSION 


The  results  of  measurements  of  bubble  resonance  f' equency 
(bubble  radius),  rectified  diffusion  rates,  and  bubble 
response  using  the  device  described  in  the  previous  chapter 
are  presented  in  this  chapter.  Before  presenting  these 
results,  a  brief  discussion  of  the  output  display  of  the 
spectrum  analyzer  is  presented.  Finally,  a  discussion  of  the 
possibilities  of  measuring  the  damping  coefficients  of  bubble 
using  the  dual  frequency  method  will  be  presented. 

A.  SPECTRUM  ANALYZER 

The  primary  detection  device  for  this  experimental  setup 
is  the  HP3585A  Spectrum  Analyzer.  This  spectrum  analyzer 
sweeps  through  a  specified  frequency  range  and  displays  the 
response  as  a  graph  of  amplitude  versus  frequency. 

With  all  transducers  off,  the  spectrum  would  be  a 
horizontal  line  with  ambient  noise  fluctuations  of  +/-  2  dBV. 
If  a  bubble  is  present  at  the  desired  position  and  both  the 
levitation  frequency  of  21.90  kHz  and  the  carrier  frequency 
of  2.65  MHz  are  being  transmitted,  the  display  under  of  the 
received  signal  under  these  conditions  is  similar  to  that 
shown  in  Figure  7.  The  spectral  lines  are  those  corresponding 
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Figure  7 .  Frequency  Spectrum  with  Bubble  Carrier 
Levitation  Transducers  on. 


to  the  frequency  summation  of  the  carrier  fl  and  the 
levitation  f3  signals  (fl+f3)  and  its  second  harmonic 
(fl  +  2*f3)  .  Figure  8  is  the  same  display  after  the  exciter 
transducer  has  been  activated.  The  spectral  line  displaying 
the  summation  of  the  carrier  frequency  and  the  bubble 
resonance  frequency  ( fl+fO )  is  easily  observed.  From 
observations,  the  bubble  resonance  frequency  is  obtained.  It 
should  be  noted  that  a  range  of  frequencies  on  either  side  of 
( fl+fO )  are  also  present  in  the  display. 

B.  BUBBLE  RADIUS 

The  frequency  of  highest  amplitude,  other  than  the 
levitation  frequency,  allows  determination  of  the  bubble's 
resonance  frequency.  From  the  resonance  frequency,  equation 
4,  and  Figure  2,  the  bubble  radius  can  be  obtained.  The 
bubble  rise  time,  described  by  equation  5,  can  be  used  to 
independently  determine  the  bubble  radius.  Figure  9  shows  the 
results  of  a  large  number  of  such  comparisons  between  the  two 
methods.  After  many  attempts  in  measuring  the  rise  time  of 
the  bubble  through  25  mm  of  water,  it  was  determined  that  time 
was  accurate  to  within  +/-5%.  This  error  is  displayed  on  the 
horizontal  axis  of  Figure  9.  Perfect  agreement  would 
correspond  to  all  the  data  points  falling  on  a  line  having  a 
slope  of  one.  The  slope,  determined  by  a  least  square  fit  to 
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Figure  8 .  Frequency  Spectrum  with  Bubble  Carrier 
Exciter  and  Levitation  Transducers  On. 
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Figure  9 .  Bubble  Radius  as  Determined  from  Two  Different 
Methods 


the  data,  is  0.99.  Therefore,  the  two  methods  agree  to  within 
approximately  1%  in  most  cases. 

C.  RECTIFIED  DIFFUSION 

In  a  second  verification,  the  dual  frequency  method  was 
used  to  measure  the  growth  of  bubbles  due  to  rectified 
diffusion.  These  measurements  are  done  for  bubbles  with 
radius  in  the  proximity  of  50  Jim  which  can  be  observed  with 
the  naked  eye.  For  such  a  bubble  of  50  Jim  radius,  equation  7 
gives  a  rectified  diffusion  threshold  of  approximately  23.0 
kPa.  Three  different  rectified  diffusion  rates  were  studied 
in  this  thesis.  The  first  is  for  the  case  of  a  negative 
rectified  diffusion  rate,  in  which  the  bubble  dissolves.  A 
pressure  amplitude  of  7 . 6  kPa  was  applied  to  the  bubble  using 
the  levitation  transducer.  Figure  10  shows  the  measured 
values  for  a  bubble  reduction  in  size  from  61.9  to  44.0  Jim  in 
a  span  of  585  seconds.  The  different  type  of  data  points 
refer  to  different  bubbles.  The  solid  line  is  the  direct 
result  of  the  integration  of  equation  4 . 

The  next  rectified  diffusion  measurement  is  for  the  case 
when  the  sound  field  is  near  the  threshold  value.  In  this 
case,  little  change  in  radius  is  expected.  Figure  11  overlays 
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21. 


the  first  five  spectrum  analyzer  displays  of  such  a  case, 
while  Figure  12  shows  the  radius  versus  time. 

The  final  rectified  diffusion  measurement  is  for  the  case 
when  the  sound  field  is  well  above  the  threshold  value. 
Figures  13  and  14  display  the  results  when  the  pressure 
amplitude  is  30.4  kPa.  Data  are  shown  for  four  different 
bubbles  in  Figure  14.  The  solid  line  in  Figure  14  is,  again, 
the  integration  of  equation  8 . 

These  results  clearly  demonstrate  the  ability  to  measure 
rectified  diffusion  using  the  dual  frequency  method.  The 
reader  should  compare  Figures  10,  12,  and  14  with  similar 
measurements  made  visually  by  Crum,  Figure  4. 

D.  DAMPING  COEFFICIENTS 

While  it  is  shown  that  the  device  is  very  effective  in 
measuring  the  properties  of  the  bubble  such  as  its  radius  and 
rectified  diffusion  threshold  and  rates,  some  very  important 
additional  information  can  be  retrieved  from  the  spectral 
response.  In  particular,  it  should  be  possible  to  retrieve 
the  frequency  response  and  thereby  determine  the  dampinc 
coefficient,  through  the  quality  factor.  Determining  a 
quality  factor  requires  a  relatively  smooth  response  which  is 
not  yet  provided  by  this  device.  The  explanation  for  this 
uneven  response  rests  in  the  fact  that  the  pressure  produced 
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Figure  12.  Bubble  Radius  Variations  at  Rectified 
Diffusion  Threshold  for  Air  Bubble  in  Water, 
f  =  21.9  kHz,  P  =  23.0  kPa. 


by  the  exciter  transducer  is  not  flat  once  it  is  incorporated 
in  the  device.  This  is  probably  the  result  of  reflection  from 
the  walls  of  the  cylinder  and  the  surface  of  the  water. 
Fortunately,  this  problem  can  be  avoided  by  using  the  small 
piezoelectric  ceramic  hydrophone  (originally  used  for 
calibration  purposes)  as  the  exciter  transducer.  Figure  15 
displays  the  difference  between  the  received  signal  amplitude 
with  the  small  ceramic  transducer  as  a  exciter  transmitter 
placed  directly  in  the  chamber  about  25  mm  from  the  position 
of  the  levitated  bubble  (heavy  line)  and  the  Panametrics  V301 
exciter  (thin  line) .  Figures  16  and  17  display  the 
improvement  accomplished  b^  such  a  change  in  exciter. 

Another  way  to  use  the  dual  frequency  method  to  measure 
the  damping  coefficient  is  to  use  equation  2  which  describes 
the  relation  between  the  acoustic  pressure  amplitude  at  the 
frequency  (fl+f2)  and  the  damping  coefficient.  Unfortunately, 
this  method  is  nonreliazable  in  this  thes:'  s  due  to  the 
inability  of  achieving  an  absolute  pressure  amplitude  reading 
at  such  high  frequency  with  the  equipment  available. 
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Figure  15.  Comparison  of  Exciters,  Small  Ceramic 
Transducers  (heavy  line) ,  Panametrics  V301 
(thin  line) . 
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Figure  17.  Bubble  Response  from  Small  Ceramic  Exciter 


V.  SUMMARY 


The  goals  of  developing  and  testing  a  new  device  for 
measuring  bubble  properties,  based  on  the  dual  frequency 
method,  have  been  accomplished.  The  experimental  setup  has 
demonstrated  the  ability  to  detect  single  bubbles  with 
resonance  frequencies  from  30  to  100  kHz.  From  these 
resonance  frequencies,  the  bubble  radius  was  computed  and  the 
accuracy  of  this  determination  was  verified  by  comparison  to 
radii  determined  from  the  rise  time  method.  This  comparison 
gave  an  average  difference  of  approximately  1%. 

The  ease  and  accuracy  in  measuring  rectified  diffusion 
rates,  negative  or  positive,  using  the  dual  frequency  method 
was  also  demonstrated. 

The  use  of  the  method  also  has  great  potential  for 
measuring  the  frequency  response  and  damping  coefficient  of 
bubbles.  The  response  in  the  sidebands  has  received  very 
little  attention  previously.  It  remains  to  be  seen  what 
information  can  be  obtained  from  responses  such  as  shown  in 
Figures  16  and  17. 

There  are  also  a  number  of  areas  for  application  of  this 
technique  such  as  to  study  the  effect  of  varying  a  fluid' s 
viscosity  or  surface  tension  on  a  bubbles  damping  constant. 
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